Magnetoelectric Behavior due to Phase Separation and Charge Locahzed States 
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We develop a theory of the magnetoelectric effects which appears in the vicinity of the second order 
phase transition in the low doped manganites as a result of the long range Coulomb forces. The long 
range Coulomb interactions is responsible for the formation of the inhomogeneous charged states 
and determines the characteristic length scales. We derive the phase diagram of the inhomogeneous 
charged states in the framework of the phenomenological theory of phase transitions. The large value 
of static dielectric function reduces characteristic value of the Coulomb energy of the inhomogeneous 
state and makes the appearance of the magnetoelectric effect possible. 
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The importance of the inhomogeneous phase segre- 
gated states in explanation of the anomalous transport 
and magnetic properties in manganites [IHll and in high 
temperature superconductors is widely discussed. 

In doped manganites the complex interactions between 
different degrees of freedom lead to unusual magnetic 
and transport properties. Moreover it was sugg ested[l3 
that spiral magnetic order observed in un-doped mangan- 
ites ReMnOa may lead to multiferroic behaviour. It was 
also discussed [ll| that the charge ordering in magnetic 
systems may cause the magnetoelectric effect. Here we 
propose that the magnetoelectric effect may appear not 
due to charge ordering but due to charge segregation, 
which appears near Coulomb frustrated second order 
phase transition. In this paper we discuss the tendency 
and conditions of the formation of inhomogeneous states 
with the spacial charge localization and the phase separa- 
tion within the phenomenological theory and clarify the 
role of the Coulomb interaction in this phenomenon. We 
demonstrate the possibility of magnetoelectric behavior 
in the low doped manganites and clarify the role of the 
Janh- Teller interaction. 

The problem of Coulomb-frustrated phase separation 
in different c harged systems is the subject of ongoing dis- 
cussion [sl-fgl. Il2|- 17|. Numerous studies have focused on 



first-order phase transitions where the charge density is 
coupled linearly to the order parameter (as an external 
field) [3, [l^ or to the square of the order parameter 
(local temperature) [l^. However, the Coulomb interac- 
tion has mostly been discussed on the qualitative level 
or as a factor of minor importance 17[, rather then an- 
alyzing the problem in general. We accurately analyze 
the total energy of the system and treat Coulomb inter- 
action exactly because it is the key factor which deter- 
mines the characteristics of inhomogeneous states. The 
idea of phase segregation and inhomogeneous charge dis- 
tribution in manganites was successfully applied for the 
description of the magnetoresistive effect in the limit of 
the percolating charged regions In this paper we 

emphasize the possibility of magnetocapacitance effect 
due to polarization of the non-percolated nano-regions in 
external electric field in the inhomogeneous state in the 



low-doped systems. 

The phenomenological approach to the theory of 
Coulomb frustrated phase transition emphasize that the 
properties the system are universal and are determined 
by the closeness to the phase transition point and by 
the dependence of the critical temperature of the phase 
transition on doping. This approach is essentially inde- 
pendent on other properties of the system. Therefore, 
the results demonstrate that the properties of the inho- 
mogeneous states are determined by the proximity to the 
phase transition and the strength of the Coulomb inter- 
action interactions but are independent on other micro- 
scopic interactions in manganites. Our estimates show 
that the Coulomb energy in the charged separated states 
is relatively small. 

We consider a doped system in the vicinity of a sec- 
ond order phase transition. We assume that the free 
carrier density p is proportional to the dopant concen- 
tration X = p/poi where po = e/a^ characteristic carrier 
density, e is elementary charge, and a is the lattice pa- 
rameter. The thermodynamic potential ^ — J 0('7, p)cPr 
describes the behavior of the order parameter 77 near the 
second order phase transition and the coupling between 
the order parameter and the charge density reads: 
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where 0o is the density of the thermodynamic potential 
in the high-temperature phase, 0,, is the density of the 
thermodynamic potential of the low-temperature phase, 
a, /?, ^ are coefficients in the expansion of the ther- 
modynamic potentials in powers of the order parameter 
{a — a {T — Tc), where Tc is critical temperature in the 
absence of doping, a — 1/C, C is the Curie constant), ^ 
is proportional to the diffusion coefficient, and H is the 
external magnetic field, (pint describes the interaction of 
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the order parameter with the charge density, ctq is the 
constant of interactions. In order to provide global sta- 
bility of the system we require ctq > 4'e describes 
the charging effects due to Coulomb interaction between 
the charge carriers, and p is the average charge density. 
We assume that Coulomb contribution is the strongest: 
1? <C a^7. The effect of average dopant concentration p 
to the thermodynamic potential in the high-temperature 
phase is included in (j)Q. The coefficient 7 is inversely 
proportional to the static dielectric constant e, 7 = Pq/s. 
As a result the coefficients in the thermodynamic poten- 
tial depend on charge density p and the term ap{po — p) 
determines the shift of the critical temperature Tc due to 
variation of the local charge density: 
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Here x — p/pQ and x = p/pq are dimensionless carrier 
densities. In the uniform case we have x = x {p = p) 
and ^u = H/2a{x,T) for T > and r?„ = -a(x,T)//3- 
H/a{x,T) for T <Tc {a{x,T) = a{T) - ax{l - x)) 

In the equilibrium state the order parameter and 
the carrier density are determined from equations: 
d(j){ri,x) / drj = 0, d4>{rj,x)/dx = pe, Me is the chemical 
potential, which can be found from the conservation of 
the total number of carriers. The second equation de- 
termines the dependence of the charge density x^ on 77 
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Substituting Eq.(3) to Eq.(l) {1} < 0^7) we obtain the 
expression of the thermodynamic potential (1) 4>ps for 
the equilibrium distribution of the charge density Xsiji): 
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The negative sign in the third term of Eq.(4) indicates 
that the uniform state may be unstable towards inhomo- 
geneous fluctuations. 

This instability leads to spatially inhomogeneous solu- 
tions. The most simple solution is spherically symmetric 
bubbles of the low symmetry phase separated by the large 
distance from each other. Thus this equilibrium inhomo- 
geneous state in the low doped samples has the spherical 
form and the distribution of charge has a shape of electric 
double- layer [3]. This solution has characteristic size i?o- 
The characteristic average charge density inside of bub- 
bles is xi and the average charge density outside of these 
regions is X2 ■ The average charge density in the system is 
X. The average value of the order parameter inside of the 
bubble rji ~ r]Q+H/2{—a{xi,T)) is relatively large, while 




FIG. 1: Magnetic field dependence (temperature dependence 
in inset) oi Ra, m and Q. 



the order parameter outside of bubbles is much smaller 
772 = H/a{x2,T) where rj^ ~ ~dt{xi,T) / p. It is clear 
that the charge is concentrated near the surface of the 
sphere. Therefore we can apply the approximation of the 
double electrical layer for evaluation of the Coulomb en- 
ergy. At large distances from the sphere the order param- 
eter and the charge density are equal to their equilibrium 
values X = X and 77 = ris{x) {rjs ~ ^72(2;) = H / a{x,T)). 
Therefore the thermodynamic potential $s of the volume 
Vo in that case has the form: 
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here $,0 = 0oK, - H^Vo{2&{x2,T))-\ A{xi) = 
TTa^{xi,T)p-^ ~ 2nH^a-\x2,T) + 0{H^), B{xi) ~ 
8TT^d-^p-\-a{xi,T) + Ht]^\xi,T) + 
H\Afi4{xi,T))-^), C(xi) ^ ^d{x ~ x^fF. 
Here we define dimensionless factor F [F — 
4 / (j)f,<Pr /^R'^d{x — xif') in order to parameterize 
the distribution of charge. F = 1/18 in the limit of the 
double electrical layer, Vq is the volume per one bubble 
(Vq = V/n where n is the number of spheres), d is the in- 
terphase boundary thickness. Here we take into account 
the strong screening of the localized charges and write 
the energy in the double electrical layer approximation. 
In that case it is proportional to di?g. If the screening 
is absent this energy is proportional to i?o 0- The 
parameters A[xi) and B{xi) (5) depend on temperature 
T and magnetic field H, and parameter C{xi) depends 
on the average charge density x. Parameters {x,T) are 
external and we find the phase diagram as a function of 
these parameters. Since x is determined by the doping, 
we describe the evolution of the properties of the system 
with doping. 

The minimum of the potential <I>s(i?o,a;i) (5) is deter- 
mined by the set of equations: d^s{Ro, xi)/dRo = and 
d^siRo,xi)/dxi — 0, which define the equilibrium val- 
ues of parameters: Rs and xi^. For the equilibrium size 
of the charged domain R^ we obtain: 



Rs = 
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FIG. 2: The temperature dependence of x ^{T) for two 
different magnetic fields. 



And xis is determined by the equation: 

(y,,Ao(xis) - A(xis)){{v{xis)A(3(xis) - A(xis) 



(y{xis) ^Vri ~ l){y{xis)yn ^ - l)B{xis)C{xis) 



= 1 



(7) 



where we use the notations: 773 — ?7o — 
H{-&{xi,T))-\ Aq{xi) = TTa^xi,T)p~\ y{xi) = 
4<j{^-2xi){xi-x){-3a{xi,T))-\ = ^rj^iSTj^)-^ 
For computer simulation we assume that X2s — 
X — X{xis — x), where A is ratio of the volume of the 
bubble (Vi) and the effective volume surrounding the 
bubble (V2), where the charge density x differs from x. 
As a result for the charge of unite bubble Q we have: 



Q{T,H)=po{x,-x)—Rl 



(8) 



Average magnetization m is determined by equation: 

47r Rl 

3 Vo 



m[T,H) ^rii — — 



(9) 



The temperature and magnetic field dependence of 
i?(r, H), Q{T, H) and m{T, H) are shown in Fig.l(a,b). 
Magnetic susceptibility x{H^T) is temperature depen- 
dent (x = dm{H) / dH) , as shown in the Fig. 2. 
Similarly we obtain the electric susceptibility: 



eeff{T,H) = kpQ{xi 



x)—Rl 
' 3 " 



(10) 



Here k is the coefficient which determines the local lat- 
tice deformation which appears due to charge localization 
and Jahn- Teller effect. Any charge displacement causes 
additional local deformation of the lattice and therefore 
protects against large charge displacements. As a result 
this local Jahn- Teller deformation may be estimated as 
djT{r) '~ kjTPo{x{r) — x) {k ^ ^jt) effective di- 
electric constant eeff{T, H) (10) is large because of large 
value of k. Importantly eeff{T, H) depends on magnetic 
field, and therefore has magnetoelectric properties. 

The condition that the inhomogeneous phase has lower 
energy then the uniform state ($s(i?0i2;i) < <&«««/) de- 
termine the region of the stability of the inhomogeneous 
phase. The equation $s(i?s,a;i) = 0o^o — /2a{x,T) 



define the upper boundary of the existence of the inho- 
mogeneous state. As a result this formula and Eqs.(6,7) 
define the upper boundary Tsi{x) in the recurrence form 
as a function of xi: 
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These equations determine the transition temperature 
Tgi to the stable inhomogeneous phase as a function of 
external parameters x. This condition means that the en- 
ergy of this inhomogeneous phase is lower then the energy 
of homogeneous state. Applying similar procedure we ob- 
tain the equation which defines the lower boundary of the 
inhomogeneous phase Ts2{x). Note that lower boundary 
Ts2 will be always close to the temperature Tcp(a;), be- 
cause there is no any gain in energy, when the bubble of 
the high-temperature phase is formed {^,^{ri = 0) = 0, 
and $,,(770) < 0). It leads to the essential difference be- 
tween the formation of the bubble of the low-temperature 
phase surrounded by the high-temperature phase and the 
bubble of the high-temperature phase surrounded by the 
low- temperature phase. The first one is energetically 
favorable and therefore the region of the existence of 
this bubbles is considerably large. It is important to 
note that Eqs.(lO) which determine the phase diagram 
does not depend on d. This fact allows us to avoid opti- 
mization of the thermodynamic potential with respect to 
d. Typical phase diagram of the inhomogeneous state is 
presented in Fig. 3. Phase transition to the nonhomoge- 
neous state represents typical first order phase transition. 
Metastable inhomogeneous phase appears at the temper- 
ature Tso which is much higher then the temperature of 
phase transition and it is shown in the phase diagram 
by dashed line. This line is determined by Eqs.(6,7) and 
A{xQf ^ AC{xo,x)B{xo,x). 

In order to make the phase separation possible it is suf- 
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FIG. 3: Phase diagram of the system with Tc = — GOiC, 
a la = 1750is:, Z = 2 10^: PP is is the region of the sta- 
ble high-temperature phase; FP is the region of the stable 
high-temperature phase with metastable bubbles of the low- 
temperature phase; PS is the region of the stable inhomoge- 
neous phase and metastable homogeneous high-temperature 
phase. 
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ficient to have only small variation of the charge density 
per unite cell in comparison with average charge density 
e{xi — x) ~ 0.1 — 0.2e. We can estimate the Coulomb 
contribution to the free energy Up as well as the energy 
gain due to formation of the low-temperature phase it,, 
(the forth and the second terms in Eq.(5)) per one unite 
cell. For the Coulomb contribution we obtain: 

Substituting e ~ 30 — 40 (the case of manganites), 
F ~ 1/18, i?o - 1 - 6mn{d - Rq), to Eq.(13) we ob- 
tain that Up is less then 0.02eV^. The energy gain due to 
formation of the low-temperature phase is u^j c± Sk^Tc 
{u,j ~ a^ia{xi,T - 0))V/3 ~ a^Tcv^{xi,T = 0)/C, 
C = ril{T = 0)/3kBNA, where Na is the Avogadro 
number). Therefore u^i — 0.03 — O.leV > Up and the 
phase separation becomes favorable. Therefore the typi- 
cal length scale of the nano-regions may be estimated as 
i?o — 1 — bnm. 

The analysis of the pair distribution function obtained 
by neutron scattering shows that the charge density in 
manganites is localized on the scale of 3 to 4 interatomic 
distances j5l.[l9|. The extra charge in that case is not more 
then 0.1 — 0.2e per unite cell. This is consistent with our 
estimates. This state is characterized as the state with 
nano-dimensional charge and phase separation. Dynam- 
ics of these charged nano-regions may lead to high value 
of the effective dielectric constant eeff{T,H) (10) in the 
low frequency range [20j . 

We have shown that the second order phase transition 
with strong dependence of critical temperature Tcp{x) (2) 
on doping is unstable with respect to the formation of the 



spatially inhomogeneous charged states. Within the phe- 
nomenological Landau theory we have shown that these 
states appears at some temperature Tsi, which is sub- 
stantially higher then the temperature Tcp (Fig. 3). As a 
result the phase transition becomes effectively first order 
phase transition. Note that the Coulomb interaction de- 
termines the charge distribution, the screening and the 
characteristic length scale of the nonhomogeneous states. 
The spatially inhomogeneous state becomes possible in 
the systems with the large dielectric constants and with 
relatively small charge density variations. We demon- 
strated that the effective dielectric constant Sef / {T, H) 
(10) is large and depends on magnetic field. 

In conclusion we underline that the localized charged 
states and the phase separation appears even in the case 
of the second order phase transition. Properties of these 
states are described within the phenomenological theory 
of the phase transitions. The Coulomb interaction deter- 
mines the spatial charge distribution, the screening and 
the characteristic length of charge localization. The in- 
homogeneous states become possible because of large di- 
electric constants and relatively small spatial variation of 
the charge density. Thise states in the low doped man- 
ganite may lead to the magnetoelectric behavior. The 
magnetocapacitance effect in the inhomogeneous state 
in low-doped systems may become the method of inves- 
tigation of the inhomogeneous charge segregated states 
including low-doped cuprates near the threshold of the 
superconductivity. 
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